Methanol-water solutions: a bi-percolating liquid mixture 
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An extensive series of neutron diffraction experiments and molecular dynamics simulations has 
shown that mixtures of methanol and water exhibit extended structures in solution despite the 
components being fully miscible in all proportions. Of particular interest is a concentration re- 
gion(methanoI mole fraction between 0.27 and 0.54) where both methanol and water appear to form 
separate, percolating networks. This is the concentration range where many transport properties 
and thermodynamic excess functions reach extremal values. The observed concentration dependence 
of several of these material properties of the solution may therefore have a structural origin. 

PACS numbers: 82.70.Uv, 83.85.Hf, 61.20.-p 



I. INTRODUCTION AND MOTIVATION 

In aqueous solutions, amphiphiles show very rich and 
interesting behaviour governed by the tendency of the 
molecules to self-organise into structures where the hy- 
drophobic regions of molecules tend to be pushed to- 
gether and away from the water, enabling the hy- 
drophilic headgroups to hydrogen-bond more easily to 
the surrounding water molecules. This results in various 
supramolecular assemblies including micelles, columnar 
phases and lamellar structures depending on concentra- 
tion and temperature. An emerging route toward the 
development, testing and refinement of detailed molecu- 
lar models of the hydrophobic interaction, hydration and 
the physics of aqueous macromolecules involves the use of 
small molecule systems (such as lower alcohols) as "pro- 
totypes" . 

Despite their structural simplicity, it is well known that 
the thermodynamic and transport properties of alcohol- 
water mixtures, such as the mean molar volume, the 
diffusion coefhcient, the compressibility and the excess 
entropy, are significantly smaller, and the viscosity sig- 
nificantly larger, than the values that might be expected 
from an ideal mixture of the pure liquids[i|,|2,yjSISIa|- 
The longstanding explanation of these effects in terms 
of an enhanced structuring of the water in the pres- 
ence of the alcohol|2| does not appear to be supported 
by modern diffraction experiments^, 0, 0, ^| and an 
alternative model is needed. Recent neutron diffraction 
studies of alcohol- water binary mixtures 's'. 'ot IT^ IT^ are 
leading to new insights into the behaviour of water near 
molecules containing both hydrophobic and hydrophilic 
groups. These have established that, in the dilute alcohol 
limit, the alcohol molecule appears to have a mildly com- 
pressive effect on the water structure, as is seen from the 
slight inwards movement of the second peak of the water- 



oxygen radial distribution function compared to the same 
function in pure water. This second peak, which occurs 
near r « 4.5 A in the OwOw radial distribution fmic- 
tion of pure water, has widely been interpreted as the 
signature of the tetrahedral ordering in water. By con- 
trast, in the opposite (concentrated alcohol) limit, the 
system segregates into what is effectively a molecular- 
scale microemulsion 8] , with methyl head groups pushed 
towards each other, and the hydrophilic hydroxy 1 groups 
forming a boundary around small pockets of a water-like 
fluid. 

These simple systems have also been the subject of 
considerable coniputational investigations. The earliest 
of these 0, 0, 0| used Monte Carlo methodologies at 
low or infinitely dilute concentrations of alcohol. De- 
spite different computational models, and some apparent 
contradictions between their results, they all found an 
enhanced cagelike structure of water around the methyl 
group, in accordance with the Frank and Evans model 7]. 
Later, MD simulations explored other mixture composi- 
tions [17L ITsI using effective potential models. Tanaka 
and Gubbins were amongst the first to highlight the 
role of the water-water interactions in discussin g aq ueous 
solutions. More recently, Meng and KoUman |2Q| have 
performed MD simulations of various solutes (including 
methanol) at infinite dilution and found that the wa- 
ter structure around the hydrophobic groups is preserved 
rather than enhanced. Laaksonen et al pll | have explored 
several concentrations though there was no specific atten- 
tion given to clustering. Ah initio simulations of alcohol- 
water mixtures have also recently been reported [23. 123l| , 
however the computational expense of these simulations 
is such that they are restricted to picosecond simulation 
on small system sizes. Nonetheless, these studies have 
also pointed to the lack of structural enhancement of the 
water surrounding the hydrophobic moiety in the alco- 
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hoi. Very recently, some of us have obtained preliminary 
results from MD simulations of an alcohol-rich methanol- 
water solution that does exhibit extreme clustering and 
micro-immiscibilitvp^ 

Given the ongoing interest in these svstemspsl |2^ . 
the availability of experimental data only at dilute and 
concentrated alcohol limits, and the apparently contra- 
dictory results from computer simulation there is a strong 
motivation to undertake a systematic survey of extended 
structure (clustering) as a function of concentration in 
the model aqueous methanol system using both experi- 
mental and simulation techniques performed at identical 
state points. We are specifically interested in exploring 
the changes in the clustering behaviour as a function of 
concentration, and the extent to which molecular dynam- 
ics simulations account for the experimental observations 
and can provide additional insight into clustering dynam- 
ics. 



II. EXPERIMENTAL METHODS 

A. Neutron diffraction experiments 

Protiated and deuteriated samples of methanol and 
water were obtained from Sigma- Aldrich and used with- 
out additional purification. Neutron diffraction mea- 
surements were performed on the SANDALS time-of- 
flight diffractometer on the ISIS pulsed neutron source 
at the Rutherford Appleton Laboratory, U.K. Samples 
were placed in flat plate cells made from a Ti-Zr alloy that 
gives negligible coherent scattering. These were mounted 
on a closed cycle refrigerator, and neutron diffraction 
measurements were made at temperatures of 293K (with 
mole fraction x = 0.05, x = 0.27, x = 0.7) and 298K 
{x = 0.54) respectively. Corrections for attenuation and 
multiple scattering were made using the ATLAS program 
suite. A further correction for inelastic scattering was 
also made[27j. The differential scattering cross-section 
for each sample was obtained by normalising to a vana- 
dium standard sample. A total of 7 samples were mea- 
sured - see Figure 1 for x — 0.54. These were respectively 
(i) CD^OD in D2O; (u) CD3OH in H2O; (iii) a 50:50 
mixture of (i) and (ii); (iv) CH3OD in D2O; (v) a 50:50 
mixture of (i) and (iv); (vi) CH^OH in H2O; and (vii) 
a 50:50 mixture of (i) and (vi). For x = 0.05 5 samples 
were measured (i),(ii),(iii),(vi) and (vii). These proce- 
dures lead to a structure factor F{Q) having the form 
FiSimiQ), SxaiQ), SxxiQ)) where Shh{Q,) gives corre- 
lations between labelled atoms and Sxh{Q,) and Sxjs:(Q) 
are the two composite partial structure factors which give 
the remaining correlations between other types of atoms 
(X) and the labelled atom type (H) in the form of a 
weighted sum of individual site-site partial structure fac- 
tors. 



B. Empirical potential structure refinement 

A mixture of methanol and water contains 6 distinct 
atomic components, namely C, O, M and H on the 
methanol molecule (here M represents the methyl hy- 
drogen atom while H represents the hydroxyl hydrogen 
atom), and Ow and Hw on the water molecule. A 
full structural characterisation of the system therefore 
requires the determination of 21 site-site radial distribu- 
tion functions, which is well beyond the possibility of any 
existing diffraction techniques by themselves. 

Therefore to build a model of the liquid structure, the 
experimental data are used to constrain a computer sim- 
ulation of the mixture. However, unlike conventional 
simulations the empirical potential used here is obtained 
directly from the diffraction data and has the effect of 
driving the structure of the three-dimensional model so- 
lution toward configurations that are consistent with the 
measured partial structure factors [28j. A total of 600 
molecules (methanol and water) are contained in a cu- 
bic box of the appropriate dimension to give the mea- 
sured density of each solution at the appropriate tem- 
perature(see Table I). Periodic boundary conditions are 
imposed. Reference interatomic potentials for water and 
methanol are taken from the literature 29? ]. A compar- 
ison between the experimentally-measured partial struc- 
ture factors and those generated from the ensemble- 
averaged EPSR configurations is shown in Fig^ 

In the present case a single set of site-site empirical 
potential coefficients was refined against the methanol- 
water data at each concentration, as well as for pure 
water [sij and for pure methanoljs^- The result is a 
set of site-site empirical potentials which are consistent 
with methanol-water solutions over the full range of con- 
centrations. Comparing the results of these simulations 
with those where the empirical potential coefficients are 
refined separately for each concentration revealed some 
discrepancies in the detail of the extracted Ow — Ow 
radial distribution function, particularly at the higher 
methanol concentrations. Clearly the diffraction data by 
themselves do not constrain this function sufficiently to 
give a completely unambiguous Ow ~ Ow radial distri- 
bution function. These discrepancies however do not af- 
fect the main conclusions of this paper, which are to do 
with the way methanol and water form distinct local clus- 
ters when mixed. We hope to present a more exhaustive 
study of the uncertainties involved in the EPSR analysis 
of molecular liquids and mixtures in a separate publica- 
tion. 



III. SIMULATION AND CLUSTER ANALYSIS 
METHODS 

We have performed classical molecular dynamics sim- 
ulations within the NVT ensemble, utilising previously 
tested intermolecular potentials for both methanol [33 
and water !34| that have been shown to predict the struc- 
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TABLE II: Parameters of the methanol- water mixtures used 
TABLE I: Parameters of the methanol-water mixtures used in the Molecular Dynamics Simulations 
in the Empirical Potential Structural Refinement 



ture and dynamics of the single component liquids well. 
Both molecules are modelled as fully flexible entities, 
with explicit potential terms for each type of atom cen- 
ter. Using the DL_POLY code ,35|, we have performed 
simulations of 2ns duration with a timestep of 0.5fs on 
systems of x=0.27, a;=0.54 and x—Q.l mole fraction. The 
parameters relevant to these simulations are given in Ta- 
ble II. All simulations were equilibrated for 0.5ns prior 
to data collection and trajectory snapshots were saved 
every O.lps for subsequent analysis. 

In analysing both the experimentally-constrained 
EPSR configurations and molecular dynamics trajecto- 
ries, an identical definition of a cluster is made based on 
bond connectivity. Specifically, water molecules are as- 
signed to the same cluster if they can be connected by a 
continuous hydrogen-bond network. The criterion used is 
that two water molecules are said to be hydrogen-bonded 
if their constituent oxygen atoms are less than RowOw 
apart, where RowOw determined as the minimum fol- 
lowing the first peak in the Ow^w P^^ir correlation func- 
tion (approximately 3.5A for both simulated and EPSR- 
fittcd functions). For the methanol molecules, clusters 
may be defined in one of two ways. When investigating 
hydrogen-bonded clusters, we use the same criterion as 
for water, i.e. if constituent methanol oxygen atoms are 
less then Roo apart. However we also investigate cluster- 
ing of methanol molecules via methyl group association. 
In this case, two methanol molecules are assigned to the 
same cluster if the C-C distance is less than the mini- 
mum following the first peak determined from the CC 
pair correlation function (which is approximately 5.7 A). 
According to this criterion, methanol molecules that are 
in contact only via their nonpolar groups are not within 
the same (hydrogen-bonded) cluster. 

When plotting cluster distribution histograms, wc plot 
the number of clusters of a size i, m{i), as a fraction of 
the total number of clusters, M, where M = J2i 'm{i) . 
Cluster lifetimes were deduced from the MD trajectories 
by calculating the average duration that a cluster persists 
with at least one participant member from the previous 
trajectory snapshot. The distinction is made between 
clusters of size 1 (i.e. with no H-bondcd neighbours) and 
those of larger sizes. 

We have investigated the dependence of cluster distri- 
butions obtained on the size of the system employed. Us- 
ing four different system sizes ranging from 700 to nearly 
18000 particles, we find that the cluster distributions and 



lifetimes are practically identical in all cases and conclude 
that our observations are not significantly dependent on 
the choice of system size. 

IV. RESULTS AND DISCUSSION 

A. Local structure 

A comparison of radial distribution fuctions g{r) for 
CC and OwOw is shown in Fig |2 A concentra- 
tion of a: = 0.27 and the pure solvents (water x = 
0.0 and methanol x = 1.0) are shown for both the 
experimentally-constrained EPSR configurations and the 
molecular dynamics simulations. Local structure there- 
fore appears similar in both the experimental and com- 
putational results. There are subtle differences between 
the experiment and simulation particularly the g{r) for 
CC which could be due to different potentials used in the 
fitting procedure and the different treatments of methyl 
group flexibility. However, it is clear that there is good 
qualitative agreement when we go on to explore extended 
structures in the solutions. 



B. Clustering and cluster lifetime 

We consider first the results of the neutron diffraction 
experiments. Visual inspection of the boxes of atoms re- 
veals significant segregation of water from methanol at 
all concentrations. An example is shown for x — 0.54 in 
Fig El Visual inspection also suggests that the methanol 
clusters do not tend to form hydrogen-bonded chains 
to the same extent as in the pure alcohol. Instead the 
methyl headgroups tend to be in contact, with the hy- 
droxyl headgroups bonding to water molecules forming 
the main boundary between methanol- and water-rich re- 
gions. This is broadly as expected for a hydrophobically- 
driven system and is what has been observed in earlier 
diffraction workp^ IT^. Similar snapshots are obtained 
from the molecular dynamics simulations(not shown). 

Further evidence of the way in which the presence of 
water affects the structure and orientation of methanol 
molecules is obtained from investigation of the size of 
hydrogen-bonded methanol clusters from the MD sim- 
ulations and radial distribution functions. Compared 
to the hydrogen-bonded network in the pure liquid, 
the cluster sizes are much smaller. For example, in 
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the methanol-rich solution (x=0.7) we find about ~ 
75 % of tlie liydrogen-bonded methanol clusters consist 
of only 1,2 or 3 methanol molecules. For comparison, 
in pure methanol the fraction of clusters of the same 
size range drops to around ~ 37 %, with the major- 
ity of hydrogen-bonded clusters containing more than 3 
methanol molecules '2^. This indicates a substantial dis- 
ruption to the methanol hydrogen-bonded network. The 
orientation of the methanol molecules is also affected; in 
the solution the methyl headgroups are pushed closer to- 
gether, as evidenced by a shift to smaller r values of the 
first peaks in gcc and qmai , where M denotes a methyl 
hydrogen, in both the EPSR procedure and MD simula- 
tions. 

The persistence of clustered structures in these systems 
is reflected in the average lifetime of clusters and single 
molecules. An estimate of these lifetimes can be obtained 
from analysis of the molecular dynamics trajectories and 
are particularly interesting in the simulation performed 
at the X = 0.7 methanol mole fraction. Firstly, we ob- 
serve that individual water molecules are short-lived and 
survive, on average, for only 2 ps before being absorbed 
into a cluster. However, in rare cases, lifetimes of 100 
ps are found. A similar result is found for the other 
cluster sizes which show average lifetimes of about 3ps 
though there are also persistent clusters surviving for up 
to 0.5 ns. The methanol hydrogen-bonded clusters, al- 
ready noted to be much reduced in size by the presence 
of the water, are extremely short-lived; most persist for 
approximately Ips with no methanol hydrogen-bonded 
structure lasting for more than 40ps. Thus the simu- 
lation results, in addition to being consistent with the 
overall structures implied by the neutron diffraction data, 
also suggest that the extended structures characterising 
the methanol-water system are very dynamic with rapid 
shedding and reforming of cluster members. 



C. Percolation 

We now explore clustering of both species quantita- 
tively as a function of concentration where an identi- 
cal definition of a cluster is made in analysing both 
the experimentally-constrained EPSR configurations and 
molecular dynamics trajectories. For water molecules the 
hydrogen-bond definition was used to designate which 
molecules belong to a given water cluster, while for 
methanol clusters the C-C distance definition was used, 
as this criterion is more indicative of the nature of the 
methanol clustering than the hydrogen bond criterion. 

The cluster size distributions as obtained from the 
EPSR ensembles (for x = 0.7, 0.54, 0.27 and 0.05 
) and molecular dynamics simulation (for x — 
0.7,0.54 and 0.27 ) are shown in Fig 0| along with the 
predicted power law UsK, s~^-^ for random percolation 
on a 3-d cubic lattice jSg. The experimental and compu- 
tational results both give similar results for the cluster 
size distribution at the three concentrations for which a 



direct comparison can be made as was implied by visual 
inspection of the EPSR and molecular dynamics struc- 
tures. 

Several "special" concentrations emerge as defined by 
changes in clustering behaviour. Specifically, x ~ 0.27 
determines the approximate alcohol concentration be- 
low which water percolates throughout the mixture while 
methanol does not - occurring instead only in small, 
isolated clusters. Above this molar fraction, however, 
methanol percolates throughout the mixture. The larger 
water clusters also percolate (at x = 0.54, see Fig0] ) 
but become increasingly isolated until they are confined 
to distinct, non-spanning clusters by x = 0.7. Thus, 
according both to the experimentally-constrained EPSR 
data and the MD simulations, in the approximate concen- 
tration range defined by 0.27 < a; < 0.54 both water and 
alcohol clusters percolate simultaneously, making this a 
bi-pcrcolating liquid mixture. 

Significantly, the mole fraction range over which si- 
multaneous two-component percolation occurs coincides 
closely with the concentration at which many thermody- 
namic properties show extrema^37, 38, 39, 40, 41, 42]. 
This suggests that the nature and extent of clustering in 
these mixtures may offer a structural explanation for the 
thermodynamic anomalies. 

In earlier work on liquids, computer simulations 
have identified percolation transitions in supercritical 
Lennard- Jones fluids lial . supercritical water j4J|, water 
in aqueous acetonitrile|43 and aqueous tetrahydrofuran 
(THE) However, to our knowledge this is the 

first report of simultaneous percolation of two fully mis- 
cible fiuids. Despite subtle variation in the cluster size 
distribution, we have obtained independent evidence of 
two-component percolation from experimental and com- 
putational methods. 



D. Dimensionality of clusters 

A further feature of the EPSR clusters is revealed in 
Fig[Slwhich shows the typical surfacWe now explore clus- 
tering of both species quantitatively as a function of con- 
centration where an identical definition of a cluster is 
made in analysing both the experimentally-constrained 
EPSR configurations and molecular dynamics trajecto- 
ries. For water molecules the hydrogen-bond definition 
was used to designate which molecules belong to a given 
water cluster, while for methanol clusters the C-C dis- 
tance definition was used, as this criterion is more in- 
dicative of the nature of the methanol clustering than 
the hydrogen bond criterion. 

The cluster size distributions as obtained from the 
EPSR ensembles (for x = 0.7, 0.54, 0.27 and 0.05 
) and molecular dynamics simulation (for x = 
0.7,0.54 and 0.27 ) are shown in Fig 0] along with the 
predicted power law Us ~ s~^'^ for random percolation 
on a 3-d cubic lattice [s^. The experimental and compu- 
tational results both give similar results for the cluster 
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size distribution at the three concentrations for which a 
direct comparison can be made as was imphed by visual 
inspection of the EPSR and molecular dynamics struc- 
tures. 

Several "special" concentrations emerge as defined by 
changes in clustering behaviour. Specifically, x « 0.27 
determines the approximate alcohol concentration be- 
low which water percolates throughout the mixture while 
methanol does not - occurring instead only in small, 
isolated clusters. Above this molar fraction, however, 
methanol percolates throughout the mixture. The larger 
water clusters also percolate (at x = 0.54, see Fig^ ) 
but become increasingly isolated until they are confined 
to distinct, non-spanning clusters hy x — 0.7. Thus, 
according both to the experimentally-constrained EPSR 
data and the MD simulations, in the approximate concen- 
tration range defined by 0.27 < x < 0.54 both water and 
alcohol clusters percolate simultaneously, making this a 
bi-pcrcolating liquid mixture. 

Significantly, the mole fraction range over which si- 
multaneous two-component percolation occurs coincides 
closely with the concentration at which man y th ermody- 
namic properties show extremajsj IssL Isol l4Ct liH l4a |. 
This suggests that the nature and extent of clustering in 
these mixtures may offer a structural explanation for the 
thermodynamic anomalies. 

In earlier work on liquids, computer simulations 
have identified percolation transitions in supercritical 
Lennard- Jones fiuids'Sl, supercritical water .44]. water 
in aqueous acetonitrile 45j and aqueous tetrahydrofuran 
(THF) pLoL l47j. However, to our knowledge this is the 
first report of simultaneous percolation of two fully mis- 
cible fluids. Despite subtle variation in the cluster size 
distribution, we have obtained independent evidence of 
two-component percolation from experimental and com- 
putational methods, e area to volume ratio (as repre- 
sented by the number of water molecules in a cluster 
which form hydrogen bonds with a methanol molecule 
divided by the number of water molecules in the clus- 
ter) for water molecules prior to the percolation tran- 
sition, together with an exactly analogous quantity for 
the methanol clusters. Clearly the ratio does not decay 
as iV"^/"^ as would be expected for a 3-dimensional ob- 
ject: the clusters appear to maximise their surface area 
by forming as many bonds as possible with methanol. 
The observed behaviour corresponds much more closely 
to a 2-dimensional object, suggesting the clusters occur 
in the form of disordered sheets or cylinders, rather than 
the sphere-like objects that might be expected in con- 
ventional micelle formation. Only at the highest water 
content, x — 0.05, do the water clusters appear to have 
adopted 3D characteristics. 

This broad conclusion is also supported by the analysis 
of the topologies of the water clusters predicted from the 
MD simulations. Figure shows the variation of the av- 
erage radius of gyration of all the clusters of size i {R{i)) 
as a function of cluster size, for three different solution 



compositions. The cluster topology can be characterised 
by a fractal dimension d, determined by a power law fit 
to the data in Figure IHI , such that R{i) (x i^"^. For 
solutions of methanol mole fractions of a;=0.7, 0.54 and 
0.3, these values of d arc determined to be 1.69, 1.89 and 
2.03, respectively. We note that as the mole fraction of 
water (and hence the proportion and size of larger clus- 
ters) increases the characteristic dimension increases. In- 
terestingly, the cluster topologies of small clusters deter- 
mined from our simulations are insensitive to changes in 
composition; all three compositions exhibit similar frac- 
tal dimension (d ~ 1.6) for clusters sizes up to i 20 . 
These results indicating the 2-d fractal dimension of wa- 
ter clusters over a range of compositions is in agreement 
with the results of a simulation study by Oleinikova et al 
j46j on the percolation of water clusters in the vicinity of 
a region of immiscibility in an aqueous solution of THF. 
Simulatancous bi-pcrcolation of both THF and water was 
found in this study, with THF percolating clusters having 
a characteristic fractal dimension of 2.5 and those of wa- 
ter approximately 1.9. Visual inspection of some of our 
percolating clusters (of both methanol and water) reveal 
that these clusters span all six faces of the simulation 
box. 



V. SUMMARY AND CONCLUSION 

Isotope-labelled neutron diffraction measurements 
analysed using the empirical potential structural refin- 
ment method have been combined and compared with 
independent molecular dynamics simulations at identical 
state points to explore structure in methanol- water mix- 
tures at several concentrations. We find that local and 
extended structures are well described by both methods 
and lead to similar conclusions. We find highly hetero- 
geneous mixing across the entire concentration range de- 
spite apparent miscibility of both components in all pro- 
portions. Extended chain, sheet and three-dimensional 
structures form depending on concentration. 

At a particular concentration regime near x — 0.27 
these structures form percolating networks for both com- 
ponents. This concentration has long been considered 
"special" as it is near the point where many transport 
coefficients and thermodynamic functions have extremal 
values. Other alcohols also show extrema of these same 
material properties (at lower mole fractions) and the 
present work suggests a structural basis for these observa- 
tions connected to the details of mixing heterogeneities. 
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FIG. 1: Typical example of the fits (lines) obtained by the 
EPSR computer simulation procedure compared to the orig- 
inal data (circles). The data shown in this case {x — 0.54) 
are the interference differential scattering cross-sections for 
the samples (i) through (vii) described under Methods. Dis- 
crepencies are observed in the low Q region. These are caused 
by difficulties in removing completely the effect of nuclear re- 
coil from the measured data. However this recoil effect is 
expected to have only a monotonic dependence on Q and so 
is unlikely to influence the model structure to any significant 
extent. 
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FIG. 2: Comparison of the radial distribution function 
g{r) as obtained from molecular dynamics(MD) time av- 
eragcs(lcft) and experimentally-constrained EPSR ensemble 
averages (right). The data shown is for the x = 0.27 mix- 
ture in relation to the pure solvents, water(x = 0.0) and 
methanol(a; = 1.0), for both the OwOw{iop) and the C-C 
radial distribution functions (bottom). 



FIG. 3: Snapshot of an experimentally-constrained EPSR 
model of the methanol- water mixture at x = 0.54 showing 

clusters of the segregated components. Methyl groups are 
shown as black spheres, large white spheres highlight the 
position of water molecules and small grey spheres denote 
methanol oxygen atoms. 




FIG. 4: Cluster size distributions for water(top) and 
methanol(below) clusters in methanol-water mixtures. For 
water molecules the hydrogen-bond definition was used to 
designate which molecules belong to a given water cluster, 
while for methanol clusters the C-C distance definition was 
used. On the left, from MD simulations with methanol mole 
fractions 0.27 , 0.54 and 0.7 and on the right from neutron 
diffraction experiments for methanol mole fractions 0.05, 0.27 
, 0.54 and 0.7 The dashed lines show the predicted cluster 
size distribution at the percolation thresholdlsj. Percolation 
in the simulated box occurs when clusters of a size close to 
the number of molecules in the simulation box form (vertical 
lines on the right hand side of the plot). 
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FIG. 5: Ratio of number of water molecules at the surface of 
a cluster (as defined by being hydrogen bonded to a methanol 
hydroxyl group) to total number of water molecules in a clus- 
ter(top). The dashed line shows the N~^^^ behaviour that 
would be expected for this ratio if the clusters grew equally 
in 3 dimensions with N the number of molecules in a cluster. 
Only for the fully percolating water cluster at x = 0.05 do 
the clusters show normal 3D behaviour. Ratio of number of 
methanol molecules at the surface of a cluster to total number 
of methanol molecules in a cluster(below). Here even at x = 
0.7 the methanol clusters do not approach full 3D behaviour. 
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Cluster size, i 



FIG. 6; Average radius of gyration of all water clusters of 
size j as a function of cluster size i, in solutions of different 
composition from MD simulations. The composition of the 
solutions, in terms of mole fractions of methanol, are x=0.7 
(squares), 0.54 (circles), 0.3 (stars). The solid and dashed 
lines represent the variation of average radius of gyration with 
cluster size of percolating clusters on 2-d square lattice (char- 
acterised by d=91/48) and 3-d cubic lattice (characterised by 
d—2.53), respectively, as determined from large lattice ran- 
dom site percolation simulations|36ll. 



